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Abstract 

43 ' 

O We consider a model for quark and lepton masses and mixings based on S4 

flavor symmetry. The model contains six Higgs doublets where three of them 
(]} ; give mass to the leptons and the other three gives mass to the quarks. Charged 

fermion and quark masses arise from renormalizable interactions while neutrino 
Majorana masses are generated through effective dimension five Weinberg op- 
erator. From the study of the minimization of the scalar potential we found a 
ly-) residual //or symmetry in the neutrino sector predicting zero reactor angle 

• and maximal atmospheric angle and for the quark sector we found a four-zero 

texture. We give a fit of the mass hierarchies and mixing angles in the quark 
sector. 
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a ■ 1 Introduction 



Quarks and leptons have very different mixing angles. A successful phenomenological 
ansatz for leptons has been proposed by Harrison, Perkins and Scott and is given by [1] 



y/2/3 l/y/3 
Utbm = ( -1/76 1/V3 -1/V2 

-l/VE i/\/3 1/V2 

which corresponds to tan 2 6* atm = 1, sin 2 #cho OZ = and tan 2 6* so i = 0.5, providing a good 
first approximation to the values indicated by current neutrino oscillation data (2j|3]. 
The third massive eigenstate is maximally mixed between // and r states and the second 
eigenstate is trimaximally mixed between e, // and r. Therefore the mixing matrix in 
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eq. ([T]) is called tri-bimaximal (TBM) mixing matrix. While the experimental mixing 
matrix for quarks is given by, see [I] 



In spite of the experimental progress so far we have no a compelling theoretical evi- 
dence regarding the flavor problem, namely why we have mixings ([1]) and (j2J) and why 
fermions masses are hierarchical. 

A possibility to solve the flavor problem is by assuming a symmetry between the 
three generations, extending the standard model with a flavor symmetry Gf. In past, 
when neutrino data was lacking, hypothesis on Gf could arise only from the quark 
sector. Successful ansatz for quarks was extended to the lepton sector, see for instance 
[5]. However recent discovery of large neutrino mixings suggest a different scenario. 
Successfully ansatz for the lepton sector can be extended to the quark sector. Tri- 
bimaximal lepton mixing can be simply derived by assuming A4 flavor symmetrjQ |6J 
and other discrete flavor symmetries give nearly tri-bimaximal mixing [HE]. Example 
of extension with v4 4 for the quark sector can be found in [9TTTT]. In this paper we 
consider the lepton and quark sectors simultaneously. There are two non Abelian 
discrete groups that are suitable for such a purpose, namely T' [12l[T3] and S 4 [H] 
since they contain singlet, doublet and triplet irreducible representations. It seems 
reasonable to consider models where quarks transform as 2 + 1 of Gf and leptons as 3 
of Gf in order to obtain large mixing in the lepton sector and small mixing between 
first and second families in the quark sector with heavy top quark mass. 

The group of permutation of four objects S4 is the minimal flavor symmetry of 
the mass matrix Mi and M v yielding TBM as shown in [15J . However was recently 
clarified in [16] that the symmetries of Mi and M v are not also symmetries of the 
Lagrangian and thus 5*4 is not special for TBM but it is simply one of many groups 
that can be used for TBM. Pioneer works using the symmetry group S4 as a family 
symmetry and deducing predictions for masses and mixings of fermions are in Ref. |17j . 
An interesting feature of 5*4 is that the neutrino mass matrix generated from a general 
dimension five Weinberg operator [18J like LL00/A where A is the cut-off scale, is 
diagonalized from TBM when is an ^-triplet = (01,02,03) that takes vev as 
(0) ~ (1, 1, 1). Differently in models with A4 , the general dimension five operator 
LL00/A is not diagonalized from TBM, see [19]. The motivation is that in A4 the 
contractions {LL)y and (LL)in break TBM, while in S4 the TBM is preserved since 
the two A4 representations 1' and 1" correspond to a one irreducible representation of 
S4, that is the doublet, and the contraction (LL) 2 preserve the TBM. 

In this paper, we study a model based on the S4 flavor symmetry, where the model 
is invariant under the Gf = S4 x Z\ x Z\ x {Z 2e x x Z 2t ) product. The quarks 
and charged lepton masses arise from renormalizable interactions while the Majorana 
neutrino masses arise from the general dimension five operator. We also study the Higgs 
potential invariant under the Gp symmetry and found the minimization conditions. 

In the next section we introduce the model, in section 3 we study the minimization 
of the potential, in sections 4 and 5 we study the phenomenological consequence of our 

A4 is the group of even permutations of four objects isomorphic to the group of symmetries of 
the tetrahedron. 
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Table 1: Quark, Lepton and scalar multiplet structure of our model, see text. 



model for leptons and quarks respectively, and in section 6 we give the conclusions. 

2 The Model 

The model consist on the flavor symmetry Gf = S4 x Z\ x Z 2 x [Z 2e x Z^ x ^2t) ; but 
in order to avoid unnecessary confusions, we will split the treatment for quarks and 
leptons. 

Consider the model defined in Table ([T|). The left-handed doublets transform as a 
triplet of S4, namely L = (L e , L^, L T ) ~ 3i and the right-handed fields e#, /xr, tr as 
singlets of S4. In the quark sector we assume the third family to transform as a singlet 
of 5*4 and the first and second families as a doublet, Qd = (Qi, Q2) an d q c D = (<1r 1 , Qr 2 )- 
We have six Standard Model Higgs doublets, = (<pi, 2 , <f> 3 ) transforming as a triplet 
under ^4, Ho = (Hi,H 2 ) and H s transforming as a doublet and a singlet respectively 
under the S4 flavor symmetry. Only the quark sector is charged with respect to the Z\ 
and Z\ symmetries. As a consequence of 5*4 x Z\ x Z\ assignment, the scalar field (f) 
interact only with charged leptons and Hd,s only with the quarks at the renormalizable 
level. 

In order to have diagonal charged lepton mass matrix we also assume extra auxiliary 
symmetries Z 2e x Z2 fJ , x Z2 T - In particular each right-handed field l c a is charged under 
the corresponding Z 2a with a = e, /i, r as well as each component of 4>. As in [20H22] 
the (Z 2 ) 3 symmetries glue each l\ with the corresponding fa. We will show below 
that (Z 2 ) 3 remove off-diagonal terms in the charged lepton sector. Since the (Z 2 ) 3 
symmetries do not commute with S4 we have to take the semidirect product^ of S4 
with (Z 2e x Z 2/1 x Z 2t ). The assignment of the charged leptons and the scalar field 
with respect to S4 x (Z 2e x Z 2/I x Z 2t ) can be summarized in the following table. 
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The 5*4 x Z\ x Z\ x (Z 2e x Z 2fl x Z 2t ) invariant Lagrangian reads 

C = L\ + £jy + C q 



2 For the use of semi-direct product in model building see for instance [22ll23| . 
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where 

Ci = y e L e e R (j>i + y^L^ R 4> 3 + y T L T T R <f) 2 ; (3) 

Cm = ^(L^)i(L^)i + ^(L 4 0) 2 (L,0) 2 + ^-(L^) 3l (L^) 3l + ^(L^) 32 (L^) 32 + 
+ ^{LiLj)x{(f)(f))i + ^(L i L,) 2 (#) 2 + ^(L^-K^k ; (4) 

£y = Y 2 d ^ D d RD H s + Y 3 d ^ s d Rs H s + Y 4 d ^d Rs H D + Y 5 d Qsq RD H D + 

+ Y 2 u ~Q d u Rd H s + Y 3 u ~Q s u Rs H s + Y^Q D u Rs H D + Y^Q s u Rd H d , (5) 

where A is a cut-off scale and cf) = ia 2 <f)* and so on. According to the S± symmetry also 
the following off-diagonal terms in the lepton sectors are allowed 

y e {L^e R (f) 3 + L T e R (p 2 ) + y^L^nfa + L T fi R <f) 2 ) + y T (L e r R (f)x + L^r R <f) 3 ), (6) 

but these terms are not invariant under the Z 2e x Z 2il x Z 2r symmetry. 

3 The scalar potential 

In our model there are 6 Higgs doublets that belong to one triplet = (fi 2 , <p 3 ), one 
doublet H D = (Hi,H 2 ) and one singlet H s representations of S , 4 . In general the Higgs 
potential can be written as 

V = V(4>) + V(H D , H.) + H D , H s ) (7) 

where V(<p) contains only ^-triplet, V(H D ,H S ) contains both S4 singlet and doublet 
scalars and VJ nt ((/>, Hp, H s ) contains only quartic terms mixing the triplet with the 
doublet and the singlet. Moreover Hp and H s transform with respect to Z 3 x Z\ and 
each component <pi of the 5*4 triplet <fi transforms with respect to Z 2e x Z 2il x Z 2r . 
Below we give the potential invariant under S4 x Z\ x Z\ and successively we consider 
the Z 2e x Z 2ii x Z 2t symmetry when we explicitly write each term in its components. 
The contributions to the Higgs Potential in eq. ([7]) invariant under the 5*4 x Z 3 x Z 2 
symmetry are given bj|f| 

V{(j>) = /i^V^ + a{[(j>^} u f + Z 2 [0%[0% + / 3 [0V] 3l [0V] 3l + 



+ ^4[WJ3 2 [WJ3 2 - 

V(H D ,H S ) = / i 2 [^iJ D ] ll + / X3[^M + a 2 ([^D]iJ 2 + a3([^D]i 2 ) 2 + 

+ a 4 [HlH D ] 2 [HlH D ] 2 + hiH^H^HlH^ + a 5 [HtH s } 2 + 

+ hiHjjHjjli^HsHs]^, 

V int (<p,H D ,H s ) = a^fl^HlH,^ +b 1 [tf(f>] ll [HlH D ] ll +b 2 [^(f>) 2 [H} ) H D ] 2 

+ c 3 [Hy} 3l [^H D } 3l +c A \Hy] 32 \^H D } 32 



3 The term proportional to le is of the form (H'H S ) Q 
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Once the symmetry Z 2e x Z 2);L x Z 2r is imposed, the contribution to the potential coming 
from Vint(0j H D , H s ) is reduced to 

M nt (0, i?£>, H a ) = [(h + c 3 - c 4 )(# 1 t #i + H\H 2 ) + ai^^] (0t 0i + t 02 + t 03); (9) 

where we can reabsorb C3 and c 4 in ftj. Similarly we reduce the other terms after 
considering the symmetry Z 2e x Z 2fJl x Z 2r , then the full scalar potential V invariant 
under S4 x Z\ x Z\ x (Z 2e x Z 2m x Z 2r ) is given by 

V = »(M\ + <h4 + <h4) + (±h + a)(M\) 2 + 

+ 2{l 2 + 5/3 - k + «)020^303 + (h + k + a)((02<4) 2 + (030 3 ) 2 ) + 

+ 2(Z 2 -Z 3 + Z4 + a)010 t l(0303>02</4) + 

+ f /^tfj + # 2 # 2 ^ 2 a 2 + + H 2 Ht) 2 a 3 + 2H x H\H 2 H\a^ 

+ [E^l + i? 2 if j) /i 2 + {H X H\ + # 2 #|) Z 5 # s #]+ 

+ ^H S H\ + a b {H s H\f + 2kH\H\H s H s + /i.e. 

+ ft!^^! + tf^X^KM + 020," + 030 3 ) + 
+ a x HlH s ((f)i4>\ + 020* + 0305) 

In the case of real vev's, that is 

(0) = (Ul, «2, U3) 

(# D > = (/l!, fc 2 ) (11) 

(H s ) = v s 

the equations of minimum are 

f£ = 2 Vl (h (h\ + h\) + h (8vl - 2 (v 2 + w 2 )) + a^ 2 + 2v\a + 2 (t; 2 + v 2 ) (l 2 + h + a)+n) 

= 0, 

g = 2t; 2 [6 1 (/2 2 + / i 2 ) + a 1 t; 2 + /i+ 

+ 2 (Z 4 (u? + uf - vf) + Z 2 (w 2 + vl) + Z 3 (-«? + v\ + 5uf ) + (w 2 + u| + uf ) a)] 

= 0, 

g = 2v 3 [b 1 (h 2 + h 2 2 ) + a 1 v 2 + f i 

+ 2 (l 2 {v\ + v 2 ) + / 4 K - v\ + w 3 2 ) + Z 3 (~vl + 5v 2 + V 2 ) + (vl + v 2 + v 2 ) a) +] 

= 0, 

(12) 
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and 

|£ = 2h x (61 (v 2 + vl + vf) + Z 5 t; 2 + 2h\{a 2 + a 3 ) + 2h\{a 2 - a 3 + a 4 ) + pi 2 ) + 2l 6 h 2 v 2 
f£ = 2h 2 (bi (vf + vl + vf) + kv 2 s + 2h\(a 2 + a 3 ) + 2h\{a 2 - a 3 + a 4 ) + A* 2 ) + 2/ 6 M s 2 
f£ = 2t; s ((h 2 + h 2 2 ) h + a 1 (v 2 + v 2 + v 2 ) + 2v 2 s a^ 3 + 2l 6 hih 2 ) = 0. 

(13) 

From the second and third equations we obtain v 2 — v 3 = v. From the first and second 
equation we found 

Vh ~ 8Z3 + 2l 4 v 3 
Vl = /, r, , r ~ = ™ ( 14 ) 

V '2 - 5/ 3 + '4 

We redefine /i 2 = and the remaining equations are rewritten as 

2v 2 (2 (l 2 - l 3 + 14 + 2/ 3 r 2 ) + (2 + r 2 ) a) + bl/i 2 (1 + <5 2 ) + /i + a lX 2 = 

bl (2 + r 2 ) v 2 + 2/i 2 (a 2 + a 3 + (a 2 - a 3 + a 4 )5 2 ) + fi 2 + (Z 5 + / 6 % 2 = 

(bl (2 + r 2 ) t; 2 + 2/i 2 (a 2 - a 3 + a 4 + (a 2 + « 3 )5 2 ) + ^2) + fe+ f )x2 = (15) 

ai (2 + r 2 ) v 2 + /i 2 (Z 5 + 2l 6 5 + l 5 5 2 ) + fi 3 + 2u^b x h\ + a lV 2 + 

+2v 2 (Z 2 (1 + r 2 ) + 2(3/3 + «) + r 2 (-/ 3 + Z 4 + a)) + /i = 

From the second and third of these equations we found two solutions for 8, one is 
5 — ±1, and the other is the one we are interested in 

hx 2 (16) 



" 2h\(2a 3 -aAy 

it is clear that in the limit Z 6 — > 0, then 5 — >■ 0. 

From the rest of equations we can determine v, v s and h\. It is straightforward 
to compute the Hessian d 2 V/duiduj, where u = (v±, v 2 , v 3 , hi, h 2 , v s ), of the Higgs 
potential. We found that for a large region of the parameter space the Hessian is 
definite positive, therefore the solution we found is a real minimum. Summarizing the 
structure of the vev's, for the doublet and the singlet of S4 we have H 1 — hi, H 2 — 5hi 
and H s — v s , so the alignment of the doublet is of the form 

H D ~(1,8). (17) 

In the limit Iq = 0, the alignment takes the form Ho ~ (1, 0). 

For the triplet of S4, Vi = rv and v 2 = v 3 = v, so the alignment is of the form 

0~(r,l,l). (18) 

Notice that in the case Z 4 = 3l 3 , implies r — 1, and the alignment takes the form 



~ 1,1,1 
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4 Lepton sector and maximal atmospheric angle 



The charged lepton mass matrix is diagonal with masses proportional to the Yukawa 
couplings y e , and y T . 

When the Higgs doublet takes vev as, see equation (1TB1) , 
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the Majorana neutrino mass matrix takes the form 
M v = h 

+ 1' 3 | 1 -l - r 2 | + u | 1 | + k | 1 o | (: 

+ k 

where U = v 2 Xi/A and r = V\/v 3 . The matrix M„ is /i <->■ r invariant therefore the 
atmospheric angle is maximal and the reactor angle is zero [24J. The solar angle depend 
from the parameter r and it is unpredicted. In the limit r — > 1 (see the potential) we 
have that x + y = w + z and the solar angle is trimaximal. 
We have three different eigenvalues: 

mi = \{™ + x + z ~ \/w 2 + x 2 + 8y 2 + z 2 + 2wz — 2x(w + z) j 
m 2 = \(w + x + z + \/w 2 + x 2 + 8y 2 + z 2 + 2wz — 2x(w + z)^ 
m 3 = z — w, 

and it is possible to reproduce the ratio a = Am 2 ol /Am 2 tm . 



(20) 



5 The quark sector 

In order to have small mixings and a hierarchical mass pattern for quarks we associate 
the quarks Ql and q c L in a 2 © li irreducible representation of 5*4. The assignment for 
particles of the model are shown in table [U From the Yukawa Lagrangian in eq. fl5]), 
once the electroweak symmetry is broken we obtain the mass matrix for the quarks 

/ Y 2 h s Y 4 h 2 \ 
M q = Y 2 h s Y^hi . (21) 
\ Y 5 h 2 Y 5 h Y 3 h s J 

From the study of the scalar potential we have the alignment 



(H D )~(1,S). 



(22) 



As we know the nearest neighbor interaction (NNI) form of the quarks matrices is in 
agreement with the experiments of quark masses and mixings, so S must be very small. 
This limit is obtained by setting / 6 = in the Higgs potential. With this fine tunning, 



the alignment for Hr> is Hp 
the form 



1, 0) and the quark mass matrices, in eq. ( 12T1) take 



M, 



u,d 



Y 2 U '% 




Y 2 U >% 



Y 3 U >% 



(23) 



Such matrices are four-zero texture and have the form of a nearest neighbor interaction, 
first proposed by Weinberg [25] and then extended by Fritzsch [26] (see also [27] and 
references therein). The mass matrices in (|23|) have factorizable phases [281129] . i.e. 



M Ujd = P^' M u ,dPR , where M has the same structure as (1231) but without phases. 
Only two combinations of phases will enter into the CKM matrix, P 



diag(l, 



pu-k pd 
L L 



* ud ). 



We rewritte the mass matrices M u d as 



M u . 



( — \ 

Vu,d 

^ b ud 

V <d vl,d I 



where 



(24) 



d 
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iPu,d + 1 -yL 
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IPv4 + 1 -Vud + R u,a 
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(25) 



and 

R*,d = v ( x + Pu, d - vtj 2 - 4 ( P u,d + < d ) + K, d yl,d 2 - 

where p u ^ and q U A defined by 



(26) 



Pa 



q 2 



ml+ral 



Pd 



m d m s 
— ^2 — 



m 2 d +m'i 



(27) 



In this case, we have 4 real parameters in each mass matrix, q Ut d, b u ^, d U A and 
y U) d- These parameters are rewritten in term of the masses of the quarks and the free 
parameters y U) d- Therefore in the mixing matrix appears, 6 masses, two real parameters 
y Ut d and the relative phases a u d = a u — a d and (3 u d = (3 U — (3d- The CKM matrix is 
then given by 

V C KM = O^PO d , (28) 
where O u d are the orthogonal matrices that diagonalize M u d via 



°ldMu,dMZd°u,d = diag(mL, m 



t,b) 



(29) 
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In this case we can fit the quark mixing angles with a very good precision, with the 
values 



y u = 0.996333 y d = 0.957981 

a ud = -2.03052 (3 ud = -1.49938 , > 

m u = 2.35634 MeV m c = 1237.37 MeV m t = 174.276 GeV 1 ' 

m d = 5.27743 MeV m s = 90.8056 MeV m b = 4243.63 MeV 



we found 



/ 0.974328 0.2251 0.00380224 
V th = 0.22497 0.973513 0.0407534 J (31) 
\ 0.00855318 0.0400268 0.999162 

which is in excellent agreement with the experimental values in eq. (J2]). The prediction 
for the Jarlskog invariant is J = 3.1 x 10~ 5 which is also in agreement withe the 
experimental central value J exp = (2.92 ± 0.15) x 10 -5 . Note that one of the phases in 
P is redundant in the sense that to make a fit of quark mixings it is enough with one 
phase, as in 



6 Conclusions 

We have studied a model based on the flavor symmetry S4 for leptons and quarks 
where all charged fermion masses arise from renormalizable Lagrangian and neutrino 
mass matrix is induced from dimension five Weinberg operator. The model contains 
six Standard Model Higgs fields transforming respectively as a triplet, a doublet and a 
singlet under S4. We have studied the minimization of the full potential. We obtain a 
H — t exchange invariant neutrino mass matrix predicting maximal atmospheric angle 
while the solar angle is undetermined. The model is compatible with normal, inverse 
and degenerate hierarchies for the neutrino masses. The quark mass matrices pattern 
have the form of a nearest neighbor interaction with four-zero texture. We give a 
numerical solution showing that it is possible to reproduce correctly the CKM mixing 
matrix within the experimental errors. It is well known that models with more than one 
Higgs SU(2) doublet may in general, have tree level flavor changing neutral currents 
(FCNC) [31 J . A complete analysis of this problem as well as a more deep analysis of 
the Higgs phenomenology, masses and decays widths, will go beyond the scope of the 
present paper, and we would like to leave this problems to a future work. An analysis 
of these problems was done in the scenario of an S3 flavor symmetry in [32] where the 
strongest constraint on the scalar masses arises from the neutral K meson mixing. 
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A The group S4 



S4 is the finite group of the permutations of four objects (for a short introduction to 
S4, see for instance [30J and references therein). S4 has 5 irreducible representations, 
two singlets li and I2, a doublet 2, and two triplets 3i and 32- 

The group £4 is defined can be defined by two generators S and T that satisfy 

S A = T 3 = (ST 2 ) 2 = 1 . (32) 

In the basis of T diagonal the generators can be written for the different representations 

as 
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In this basis the rules for the non trivial products of two irreducible representations 
are: 
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The multiplication rules with the 3- dimensional representations are the following: 

li ~ ax/3i + a 2 3 + a 3 2 



a 2 02 + ai fa + ai 3 0i 

01303 + aii p 2 + Ol 2 0l 



3i © 3i = 3 2 © 3 2 = li © 2 © 3i © 3 2 



with 



3i © 3 2 



3i 



with 



2«i/3i - a 2 3 - a 3 2 
3 1 ~ } 2a 3 /3 3 - ai/3 2 - a 2 0i 
2a 2 2 - asi0 3 - a 3 0i 

O120Z ~ OL 3 2 
01102 - 01201 
OL 3 0i - Oll0 3 



oii0i + a 2 3 + a 3 2 

01202 + 01103 + Ol 3 01 
-a 3 3 - 01x02 - 01201 



3i 




2oii0i - a 2 3 - a 3 2 
2a 3 3 - ai0 2 - ol 2 0i 

2a 2 02 - Oli03 - 01301 

The products of two irreducible representations Ax B is different from the Ax B since 
the T and the S generators are complex. When we multiply the complex conjugate 
of a two-dimensional representation with another two-dimensional representation, we 
have to interchange the indices 1 -H- 2 of the complex conjugate representation, that 
is, for instance in 2 © 2 we interchange ai -H- a 2 and the product is given by 



2 with < 



11 ~ a 2 2 + ai0i 

1 2 ~ a 2 2 - oii0i 



01102 

a 2 0i 

When we multiply the complex conjugate of a three-dimensional representation 
with another representation, we have to interchange the indices 2 -B- 3 of the complex 
conjugate representation, that is, for instance in 3©3 we interchange a 2 -H- a 3 . Similarly 
in the case of 2 © 3 we have to interchange the indices 1 -H- 2 for the doublet, that is 
oil oi2- 
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